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ABSTRACT 

Time-distance analysis and acoustic imaging are two related techniques to probe the 
local properties of solar interior. In this study, we discuss the relation of phase time and 
envelope time between the two techniques. The location of the envelope peak of the cross 
correlation function in time-distance analysis is identified as the travel time of the wave 
packet formed by modes with the same u )/L The phase time of the cross correlation 
function provides information of the phase change accumulated along the wave path, 
including the phase change at the boundaries of the mode cavity. The acoustic signals 
constructed with the technique of acoustic imaging contain both phase and intensity 
information. The phase of constructed signals can be studied by computing the cross 
correlation function between time series constructed with ingoing and outgoing waves. In 
this study, we use the data taken with the Taiwan Oscillation Network (TON) instrument 
and the Michelson Doppler Imager (MDI) instrument. The analysis is carried out for 
the quiet Sun. We use the relation of envelope time versus distance measured in time- 
distance analyses to construct the acoustic signals in acoustic imaging analyses. The 
phase time of the cross correlation function of constructed ingoing and outgoing time 
series is twice the difference between the phase time and envelope time in time-distance 
analyses as predicted. The envelope peak of the cross correlation function between 
constructed ingoing and outgoing time series is located at zero time as predicted for 
results of one-bounce at 3 mHz for all four data sets and two-bounce at 3 mHz for two 
TON data sets. But it is different from zero for other cases. The cause of the deviation 
of the envelope peak from zero is not known. 
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Soares (Instituto Astrofisica de Canarias, Observatorio del Teide, Tenerife, Spain); Guoxiang Ai and Gwo-Ping Wang (Huairou Solar 
Observing Station, Beijing Observatory, Beijing, China); Philip Goode and William Marquette (Big Bear Solar Observatory, New Jersey 
Institute of Technology, Newark, NJ 07102, U.S.A.); Shuhrat Ehgamberdiev and Oleg Ladenkov (Ulugh Beg Astronomical Institute, 
Tashkent, Uzbekistan). 
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1. Introduction 

Time-distance analysis and acoustic imaging are 
two related techniques to probe the local properties 
of the solar interior. A resonant p-mode is trapped 
and multiply-reflected in a cavity between the surface 
and a layer in the solar interior. The acoustic sig- 
nal emanating from a point at the surface propagates 
downward to the bottom of the cavity and back to 
the surface at a different horizontal distance from the 
original point. Different p-modes have different paths 
and arrive at the surface with different times and dif- 
ferent distances from that point. The modes with the 
same angular phase speed uj/1 have approximately the 
same ray path, where u> is the mode frequency and l is 
the spherical harmonic degree. In time-distance anal- 
ysis, the travel time of a wave packet formed by the 
modes with the same uj/1 versus the travel distance is 
determined from the cross-correlation function of sig- 
nals measured at two locations (Duvall et al. 1993). 
The phase time of the cross-correlation can also be 
determined (Duvall et al. 1996). The perturbation in 
phase time due to local inhomogeneities can be mea- 
sured and inverted to infer the flow, magnetic field, 
and temperature below the solar surface (Jefferies et 
al. 1994, D’Silvaet al. 1996, Duvall et al. 1996, Koso- 
vichev 1996, Braun 1997, Bogdan et al. 1997, Koso- 
vichev & Duvall 1997, more updated references?). 

In acoustic imaging, the signals measured at the 
solar surface are coherently added, based on the time- 
distance relation, to construct the signal at a point on 
the solar surface or in the solar interior (Chang et al. 
1997, Lindsey &; Braun 1997, Chen et al. 1998, Braun 
et al. 1998, Chou et al. 1999). For a target point on 
the surface, one can use the measured time-distance 
relation. For a target point below the surface, one 
has to use the time-distance relation computed with 
a solar standard model. The constructed signal con- 
tains information on intensity and phase. The phase 
of a constructed time series contains information on 
the phase change as the wave packet propagates from 
the target point to the observed point. This phase 
change relates to the phase time and envelope time 
in the time-distance analysis. In this paper, we study 
the relation of phase time and envelope time between 
time-distance analysis and acoustic imaging analysis. 

The phase of the constructed time series in acous- 
tic imaging depends on the time-distance relation 
used. If the travel time used to construct the sig- 
nals is shifted by a constant, the constructed time 


series is also shifted by the same amount (Chen et 
al. 1998, Chou et al. 1999). In this study, we use the 
measured travel time (envelope time) in time-distance 
analysis to construct signals in acoustic imaging. The 
measured time-distance relation is sensitive to the fre- 
quency (Jefferies et al. 1994, D’Silvaet al. 1999). To 
eliminate the frequency dependence, we use the same 
data set to measure the time-distance relation and 
to construct signals at the solar surface. Also, when 
comparing MDI and TON data, we isolate the same 
frequency bands for analysis. 

2. Observations and Data 

We use the helioseismic data taken with the TON 
instruments and the SOI/MDI instruments aboard 
SOHO. The TON data are full-disk K-line images 
recorded by a 16-bit 1080 x 1080 water-cooled CCD. 
The images are taken at a rate of one image per 
minute. The diameter of the Sun is about 1000 pix- 
els. The measured amplitude of intensity oscillations 
is about 2.5%. A discussion of the TON project and 
its instruments is given by Chou et al. (1995). The 
preliminary data reduction of the TON data, such as 
flatfielding and registration, is given by Chen, Chou, 
& TON Team (1996). Two 512-minute time series on 
1996 August 3 are used in this study. One is taken 
in Tenerife, defined as TON-1; another is taken in 
Tashkent, defined as TON-2. The overlap of the two 
time series is 172 minutes. 

The SOI/MDI experiment generates Dopplergrams, 
magnetograms and continuum images at different 
scales. A summary of the SOI/MDI project was given 
by Scherrer (1995). The MDI data used in this study 
are 1024 x 1024 full-disk Doppler images with one 
minute cadence. We use two 1024-minute time series 
on July 6 of 1996 (defined as MDI-1) and July 11 of 
1996 (defined as MDI-2). (Dean: I only used the first 
512 minutes for the time-distance analysis) 

Each observed full-disk image is transformed into 
sin 9-<f) coordinates with the GRASP package devel- 
oped by GONG, where 0 and <j) are the latitude and 
the longitude, respectively, in a spherical coordinate 
system aligned along the solar rotation axis. The so- 
lar differential rotation at the surface is removed with 
an observed surface differential rotation velocity (Lib- 
brecht & Morrow 1991). The data is filtered with a 
Gaussian filter of FWHM=lmHz centered at 3 and 
4 mHz. Details of data preparation are discussed in 
Chou et al. (1999). 
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3. Phase Time and Envelope Time in Time- 
distance Analyses 

The cross-correlation function C(r, A) between a 
central point and an annular region at an angular 
distance of A is computed with 

C(r,A)=/ f(t) ■ f(t + t, A) dt (1) 
Jo 

where f(t) is the signal measured at the central point, 
and f(t -j- r, A) is the signal averaged over an annular 
region at an angular distance of A from the central 
point. The width of the annular region is 0.35° for 
the TON data and 0.12° for the MDI data. This 
procedure is repeated for different central points in a 
region of 87° x 75° at the disk center. The average 
cross-correlation function C(r, A) is obtained by aver- 
aging over this region. For a fixed angular distance A, 
C(r, A) has the approximate form of a Gabor wavelet, 

G(r) = A cos(27rj/[r - r pft ]) • exp(- ^— ) , (2) 

where A, a and v are the amplitude, width and fre- 
quency respectively. The location of the peak of 
the envelope, r cn (envelope time), is identified as the 
travel time of the wave packet. The phase time, r p h> 
is undetermined to multiples of the period, i/” 1 . The 
measured average cross-correlation function at a fixed 
A, C(t } A), is fitted to the Gabor wavelet in equation 
(2) to obtain the envelope time r en and the phase time 
r p h . This procedure is repeated for different A to ob- 
tain the time-distance relations, r en versus A and r p 
versus A. 

The above analysis is carried out for two TON 
data sets (TON-1 and TON-2) and one MDI data set 
(MDI-1). The range of A is 1.75 - 24.85° for TON 
data and 1.08 — 24.96° for MDI data. The measured 
time-distance relations (both r cn and r p h) from the 
three data sets are very close. Examples of measured 
time-distance relations from the TON-1 and MDI-1 
data sets at 3 mHz are shown in Figures 1 and 2, re- 
spectively. The upper solid curve is the phase time 
r p h\ and the lower solid curve is the envelope time 
r en . Usually T p h can be determined more accurately 
than r en ■ The measured r p h versus angular distance 
is rather smooth and can be fitted well with a fifth- 
order polynomial. The fluctuation of measured r en 
is large, and it is difficult to be sure that the fit is 
good. In this study we will use the relation of r en 
versus angular distance to construct acoustic signals 


in acoustic imaging. Thus we have to find a better 
way to measure r e „. 

The difference between r ph and r en increases slowly 
with angular distance. The increase is approximately 
linear. An example (TON-1 at 3 mHz) is shown in 
Figure 3. The large fluctuation is caused by the inac- 
curacy in determining r en . The dashed line is a linear 
fit. The results from the other two data sets also have 
this property: r p / » — r en increase slowly and linearly 
with angular distance. Since measured r p /, — r en can 
be fitted well by a linear function, we subtract this lin- 
ear function from the fitted r ph to obtain r en , which 
corresponds to the dashed lines in Figures 1 and 2. It 
is noted that the value of r p ^ — r en is strongly depen- 
dent on the frequency. But it always increases slowly 
and approximately linearly with angular distance. 

The modes with the same w/l have approximately 
the same ray path, and form a wave packet. If the 
power distribution of modes is sharply peaked around 
some value wq and the second derivative of k with re- 
spect to u> is small, with the expression for wavefunc- 
tions in D’Silva (1996) and the technique in Jackson 
(1975), the wavefunction of the wave packet can be 
expressed as 


*(r,0 


i(u»oT(ro,r)— f kods—a) 


*(ro,<- T(? 0 j f)) 


( 3 ) 


where T(r 0 , r) = ff o {dk/duj)\o ds is the travel time of 
the wave packet from ro to r along the ray path. The 
term J ^ k^ ds is the phase change accumulated along 
the wave path from ro to r. The additional phase a 
accounts for the phase change at the boundaries of 
mode cavity. Equation (3) shows that, apart from an 
overall phase factor, the wave packet approximately 
maintains its form. 

The cross-correlation function of signals at r 0 and 
r is computed as equation (1) 


C(r,r 0 ,r) 

J *(r 0 ,t) ■V'(f,t + r)dt 


- i(u>oT— f kods — a) f 

J '° 'I *(r 0 ,t) **{r 0 ,t+T-T)dt 

-i(u>oT— f k 0 ds — a) ~ / 

J '° C(r — T, ro) (4) 


where C(r — T, ro) is the auto-correlation function of 
the signal at r 0 , which is peaked at r = T. Thus 
the envelope of C(r, r 0 ,r) is peaked at r = T. The 
difference between envelope time and phase time of 
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C(r, r 0 ,r) is 



where v p and v g are the unit vectors of phase veloc- 
ity and group velocity, respectively. The difference 
between the phase velocity and the group velocity is 
small except near the solar surface. 

4. Phase Time and Envelope Time in Acous- 
tic Imaging 

The acoustic signals at a point on the solar sur- 
face or in the solar interior can be constructed with 
the technique of acoustic imaging (Chang et al. 1997, 
Chou et al. 1999). Since the goal of this study is 
to compare phase information in time-distance anal- 
ysis and acoustic imaging, we will use the measured 
time-distance relation to construct acoustic signals at 
the solar surface. To study the phase information 
of constructed time series, we compare the time se- 
ries constructed with outgoing waves and the time 
series constructed with ingoing waves by computing 
the cross correlation function of two time series (Chen 
et al. 1998). The outgoing time series is constructed 
with the ordinary time-distance curve; and the ingo- 
ing time series is constructed with the time-reversed 
time-distance curve. Since our data is filtered with a 
Gaussian filter, the cross correlation function has ap- 
proximately the form of a Gabor wavelet as in equa- 
tion (2). The phase time T ph and envelope time r' n of 
the cross correlation function can be be determined 
by fitting it with the Gabor wavelet. Here we use the 
superscript prime to distinguish from the phase time 
and envelope time in the time-distance analysis. The 
phase time is also undetermined to multiples of 
the modulation period. 

In this study, we use the envelope time measured in 
time-distance analysis to reconstruct acoustic signals. 
The envelope time measured in time-distance analysis 
is interpreted as the travel time of a wave packet. 
Thus T ~ r en in equation (3). To construct the signal 
at ro and t = 0, we collect signals measured at r and 
t = T, which is ^(r, T), for various r and T, where T 
and |r — ro| satisfies the time-distance relation. ^From 
equation (3), we have 

. _ i(u>c>T(ro ,r) — f * ko da — a) 

4f(r,T) « e Jr o ^(ro,0) 


= e -*^o( T ph -r. n )^ r ( fTQ 0) (6) 

Equation (6) shows that the collected signals dif- 
fer from the original signal only by a phase factor 
exp(-iuo(T p h — r eri )) This phase varies slightly with 
distance. If the range of distance in the coherent sum 
is not large, the variation of this phase with distance 
is negligible. 

Based on the above discussion, we have two pre- 
dictions about the cross-correlation function of con- 
structed ingoing and outgoing time series, if we use 
the envelope time r en measured in time-distance anal- 
ysis to construct signals in acoustic imaging. First, 
the envelope time of the cross-correlation r ( en is zero. 
Second, the phase time of the cross-correlation Tp h is 
twice Tp/, — r en . This is because the phases of con- 
structed outgoing time series and ingoing time series 
are different from the phase of the original time series 
by T p h — r en and — (r p /, — r en ), respectively. We will 
discuss the comparison of these two predictions with 
observations. 

We have analyzed two TON data sets (TON-1 and 
TON-2). For each data set, we use the time-distance 
relation measured from this data set to construct the 
signals. We have also analyzed two MDI data sets 
(MDI-1 and MDI-2). But we use the time-distance re- 
lation measured from MDI-1 to construct the signals 
for MDI-1 and MDI-2. The range of time-distance 
curves used to construct signals is 3.4 — 7.5° for TON 
data, and 1.8 — 5.5° for MDI data. The target re- 
gion is a quiet region near the disk center. Its size is 
17.6 x 11.7° for TON data and 12.0 x 12.0° for MDI 
data. The regions used to construct signals are also 
in the quiet Sun. 

In this study, we use two different ways to deter- 
mine T f en and r ph- In the first way, we determine 
the parameters of the cross-correlation function of 
constructed ingoing and outgoing time series at each 
point by fitting it to a Gabor wavelet, and then av- 
erage the parameters over the target region. In the 
second way, we average the cross-correlation function 
over the target region, and then determine the param- 
eters by fitting the average cross-correlation to a Ga- 
bor wavelet. The two ways yield very close results. In 
the tables, we give only the result from the first way. 
The parameters of the cross-correlation function from 
the four different data sets are shown in Tables 1 and 
2, which correspond to 3mHz and 4mHz, respectively. 
The third column of the tables is the amplitude A of 
the envelope of the cross correlation function. The 
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fourth column is the modulation period, i/~ l . The 
fifth and sixth columns are the envelope time r^ n and 
the phase time T ph of the cross correlation function, 
respectively. The last column is the predicted phase 
time, which equals to 2N s (r p h - r en ), where N s is the 
number of the skip (bounce). The value of r pi i, — r en 
quoted here is the average over the range of angular 
distance used to construct signals. 

The acoustic signals can be constructed with the 
signals after more than one bounce from the target 
point. The time-distance curve of multiple bounces 
used in this study is obtained by multiplying both 
measured travel time and distance of one-bounce by 
the number of bounces. The results of 2-bounce and 
3- bounce are also shown in tables 1 and 2. 

The value of r^ n is close to zero for 1-bounce results 
of all four data sets and 2-bounce results of TON-1 
and TON-2 at 3 mHz as predicted. But it is not 
close to zero for the other cases. We will discuss the 
possible causes for the large r' n in these cases. On the 
other hand, the measured r ph is close to the predicted 
value for all cases we have studied, except one case 
(2-bounce, TON-2 at 3 mHz). We conclude that the 
measured r ph is in good agreement with the predicted 
value. 

From equation (3), the constructed time series de- 
pends on the time-distance curve used in acoustic 
imaging. If the travel time used is shifted by a con- 
stant A i, the constructed outgoing time series will 
shift by A tf, and the constructed ingoing time series 
by — At. The cross correlation function of ingoing 
and outgoing time series will shift by — 2A t. It leads 
to that both phase time and envelope of the cross 
correlation function shift by — 2A t. Therefore, it is 
unlikely that the deviation of measured r' n from the 
predicted value is caused by the inaccuracy of mea- 
sured envelope time r en used in acoustic imaging, be- 
cause measured r ph is rather close to the predicted 
value. 

The determination of r' n in the fit is less accurate 
than T ph and it is difficult to estimate the error of 
r' n . But it is unlikely that the error in the fit is large 
as the value of measured r' n , such as 2 minutes, in 
some cases, though we can not totally rule out this 
possibility. 

The power spectra of filtered data are not exactly 
a Gaussian, though we use a Gaussian filter (D’Silva 
et al. 1999). The power spectra of constructed sig- 
nals have a further distortion from a Gaussian because 


different frequencies have different lifetimes. Thus 
the envelope of the cross-correlation function of con- 
structed outgoing and ingoing time series is not ex- 
actly a Gaussian. Fitting the cross-correlation func- 
tion to a Gabor wavelet may cause some error in de- 
termining r' n . But we believe this error can not ac- 
count for the large value of measured r' n alone. 

5. Discussion 

This study confirms that the phase time r ph of the 
cross-correlation function of constructed ingoing and 
outgoing time series in acoustic imaging is equal to 
twice the difference between phase time r p h and en- 
velope time r en in time-distance analyses. The dif- 
ference between phase time and envelope time con- 
tains two information: 1) the difference between the 
phase velocity and the group velocity and 2) the phase 
change at the boundaries of mode cavity. 

If the target point is located inside a sunspot, the 
presence of magnetic field would modify both the 
phase time and the envelope time. First, the pres- 
ence of magnetic field would modify the dispersion 
relation of waves and would change the difference be- 
tween the phase velocity and the group velocity. Sec- 
ond, the magnetic field would modify the depth of the 
upper boundary of mode cavity. It makes both enve- 
lope time and phase time change. Third, the mag- 
netic field would also modify the physical conditions 
at the upper boundary of mode cavity. It would mod- 
ify of the phase change as the wave packet reflected 
at the boundary. Fourth, the phase of a wave packet 
would change as it goes across lateral boundaries of 
the sunspot. Fifth, the flow caused by the presense 
of the sunspot would change both enevlope time and 
phase time. But the changes in envelope time and 
phase time are equal. Therefore, the comparison of 
deviations of measured T ph and r' n inside a sunspot 
from the values in the quiet Sun provides information 
on the structure of the sunspot. 

We are grateful to the GONG Data Team for pro- 
viding the software package GRASP. DYC and the 
TON project were supported by NSC of ROC un- 
der grants NSC-87-21 12- M-007-044. TLD and the 
MDI project were supported by the SOI-MDI con- 
tract NAG5-3077 at Stanford University. 
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Table 1 

Phase time and envelope time in acoustic imaging (3mHz) 


data 

no. of skip 

amplitude 

period (min) 

T* a 
'en 

T t b 
T P h 

predicted r^ c 

TON-1 

1 

0.55 

5.27 

-0.05 

0.75 

0.80 

TON-1 

2 

0.57 

5.29 

0.05 

1.62 

1.60 

TON-1 

3 

0.56 

5.33 

0.14 

2.45 

2.41 

TON-2 

1 

0.51 

5.27 

-0.07 

0.77 

0.82 

TON-2 

2 

0.55 

5.27 

0.02 

1.67 

1.64 

TON-2 

3 

0.55 

5.31 

0.11 

2.53 

2.46 

MDI-1 

1 

0.72 

5.45 

0.01 

0.81 

0.83 

MDI-1 

2 

0.67 

5.51 

0.47 

1.69 

1.67 

MDI-1 

3 

0.57 

5.58 

0.22 

2.53 

2.50 

MDI-2 

1 

0.73 

5.46 

0.03 

0.81 

0.83 

MDI-2 

2 

0.68 

5.53 

0.48 

1.69 

1.67 

MDI-2 

3 

0.59 

5.60 

0.23 

2.54 

2.50 


“envelope time, in units of the period 
b phase time, in units of the period 
C predicted phase time, in units of the period 
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Table 2 

Phase time and envelope time in acoustic imaging (4mHz). 


data 

no. of skip 

amplitude 

period (min) 

r' a 

en 

r' b 
T P h 

predicted r ph c 

TON-1 

1 

0.46 

4.50 

0.22 

1.52 

1.53 

TON-1 

2 

0.40 

4.55 

0.49 

3.19 

3.07 

TON-1 

3 

0.41 

4.50 

0.30 

4.67 

4.60 

TON-2 

1 

0.42 

4.52 

0.21 

1.47 

1.49 

TON-2 

2 

0.41 

4.57 

0.50 

2.35 

2.98 

TON-2 

3 

0.42 

4.54 

0.36 

4.57 

4.48 

MDI-1 

1 

0.64 

4.46 

0.14 

1.36 

1.36 

MDI-1 

2 

0.49 

4.56 

0.28 

2.76 

2.73 

MDI-1 

3 

0.38 

4.55 

0.20 

4.18 

4.09 

MDI-2 

1 

0.64 

4.46 

0.15 

1.36 

1.36 

MDI-2 

2 

0.49 

4.56 

0.33 

2.77 

2.73 

MDI-2 

3 

0.39 

4.54 

0.28 

4.19 

4.09 


Envelope time, in unit of period 
b phase time, in unit of period 
predicted phase time, in unit of period 






Fig. 1. — Phase time and envelope time of cross- 
correlation function versus angular distance measured 
from the TON-1 data set at 3 mHz. The upper solid 
curve is the phase time r P h , and the lower solid curve 
is the envelope time r en . The dashed curve is a fit to 
the measured envelope time with a method described 
in the text. 


Fig. 3. Difference between phase time and enve- 
lope time versus angular distance measured from the 
TON-1 data set at 3 mHz. The dashed line is a linear 
fit. 


Fig. 2. — Same as Figure 1, but measured from the 
MDI-1 data set. 
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